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Abstract 

Absorption of infrared and visible radiation from stellar emission spectra by the negatively 
charged hydrogen ions H“ is considered. The explicit formula for the photodetachment cross- 
section of the negatively charged hydrogen ion(s) is derived. Photodetachemnt cross-sections of 
the °°H“, (or T“), (or D“) and ions are determined to high accuracy and for a large 
number of photo-electron momenta/energies. We introduce criteria which can be used to evaluate 
the overall quality of highly accurate wave functions of the hydrogen ion(s). One of these criteria 
is based on highly accurate calculations of the lowest order QED corrections in the negatively 
charged hydrogen ions, including (protium), (deuterium), (tritium) and model ion 
with the infinitely heavy nucleus An effective approach has been developed to calculate 

three-body integrals with the Bessel functions of different orders. Some preliminary evaluations 
of the phototdetachment cross-sections of the negatively charged hydrogen ions are performed. 
Inverse bremsstrahlung in the field of the neutral hydrogen atom is briefly discussed. 
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I. INTRODUCTION 


As is well known (see, e.g., fl-B) the negatively charged hydrogen ions (protinm) 
and (deuterium or D“) are of great interest in Stellar Astrophysics. Indeed, the nega¬ 
tively charged hydrogen ions substantially determine the absorption of infrared and visible 
radiation in photospheres of all stars with temperatures bounded between Tmax ~ 8,250 
K (late A-stars) and T^in ~ 2,750 K (early M-stars). The main contribution to the light 
absorption by the negatively charged hydrogen ion H“ comes from its photodetachment 

H- + hz/ = H + e" (1) 

where the notation hv designates the incident light quantum, while e~ means the free elec¬ 
tron, or photo-electron, for short. The notation h stands for the Planck constant. Note that 
each of the negatively charged hydrogen ions, i.e. the (or D“), (or T”) and 


has only one bound (ground) 


of this two-electron ion (see, e.g., jsj). The ground state is a singlet state, i.e., the total 
electron spin S in this state equals zero, and therefore, 25-1-1 = 1. Here and everywhere 
below in this study the notation stands for the model atom/ion with the inhnitely 

heavy nucleus. The total non-relativistic energy of the ground 1^5—state in the ion 
is -0.527751 016544 377196 59055(5) a.u. j^, while the total energy of the hydrogen atom 
with the inhnitely heavy nucleus (°°H) is -0.5 a.u. Therefore, the binding energy of the 
ion is Xi ~ 0.027 751 a.u. ~ 0.755 143 9 eV and the corresponding ionization frequency 
z/i = ^ is located in the far infrared part of the spectrum. Note that the energies of ‘visible’ 
light quanta are 2 - 3 eV, while the total energy of the ground state of hydrogen atom 
is ~ 13.605 eV. Light quanta with such energies (~ 13.605 eV) correspond to the vacuum 
ultraviolet region. 

The hnal atomic state which arises after photodetachment of the negatively charged 
hydrogen ion, Eq.([T]), includes the neutral hydrogen atom H and ‘free’ electron e~. The 
‘free’ electron moves in the held of the central hydrogen atom. The kinetic energy of this 
electron depends upon the total energy of the incident photon and can, in principle, be 
arbitrary. Formally, all radiation quanta with the frequencies i' > I'l = ^ can produce 
photodetachment of the negatively charged hydrogen ion H~. However, it can be shown 
(e.g., see below) that the photodetachment cross-section rapidly (exponentially) vanishes, if 
the photon energy substantially exceeds numerical values of the ionization potential of the 


^5—state with L = 0, where L is the angular momentum 
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H“ ion. Another reason which allows one to ignore the role of high-energy photons in the 
H“ ion photodetachment follows from stellar Astrophyscs. Indeed, if we assume that the 
actual stellar photosphere is at thermal equilibrium with T ^ 4000 - 8000 K, then the part of 
‘high-energy’ photons which can produce excitations of the hnal hydrogen atom is extremely 
small. For instance, to produce photodetachment of the H“ ion with the formation of the 
hnal H atom in the 2s—state the ‘photon temperature’ Tf must be > 78944 K. Analogous 
‘photon temperature’ to form the hnal H atom in the 3s—state is around 118,416 K. For 
our Sun the fraction of the ‘hot’ thermal photons with T > 78944 K is very small (< 1 ■ 10“^ 
%). Even in photospheres of late A—stars the fraction of ‘hot’ thermal photons is less than 
0.5 %. In reality, such photons become noticeable in photospheres of B-stars with surface 
temperature T > 20, 00077. However, in classical B-stars the photodetachment of the 
ions does not play any role and most of the emitted radiation is absorbed by neutral hydrogen 
and helium atoms. This explains why the photodetachment of the H“ ion in our Sun and 
in the late A-stars and in all F, G and K stars and in some other stars always proceeds 
with the formation of the hnal H atom in its ground (bound) Is—state. In the hot O- and 
B-stars, with surface temperatures Tg > 28,500 K and Tg > 11,000 77, respectively, the 
absorption of radiation quanta is mainly related to the transitions between different bound 
states in helium and hydrogen atoms. The negatively charged H“ ion does not exist at such 
temperatures as a stable (or bound) system. 

This situation can be different in the Be-stars, i.e. in the B-stars with the rapidly rotating, 
low-dense disk of hydrogen atoms [l|, j^. Illumination of this disk by radiation emitted by 
the central hot star produces excitation of hydrogen atoms in the stellar disk [^. Excitation 
of hydrogen atoms and following emission of absorbed radiation generates very special optical 
spectra of the Be-stars, which contain thermal (continous) spectra from the central star (or 
central cluster) and linear spectra which corresponds to the transition between excited states 
in the hydrogen atoms from the disk. For such stars photodetachment of the H“ ions with 
the formation of the hnal hydrogen atoms in their excited states may be important in the 
outer areas of stellar disk. However, for all Be-stars photodetachment of the H“ ions may 
proceed in the spatial areas located at the large distances (> 150,000,00 km) from the hot 
central star, i.e. in areas which are far away from stellar photospheres. Below in this study 
we shall not discuss photodetachment of the H“ ion at such ‘exotic’ conditions. 

On the other hand, let us note that the H“ ions play practically no role for the cold 
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M and N stars with Tg < 2,700 K where the absorption of radiation by various molecular 
species and atoms of metals becomes important. However, in the late F, G and early K stars 
the absorption of infrared and visible radiation by the negatively charged hydrogen ions is 
maximal. This includes our Sun which is a star of spectral type G2. The great role of the 
H“ ions for our Sun was suggested by R. Wild in 1939 (see discussions and references in 

n n n 

l|, [2[ and [^). An effective absorbtion of large amount of infrared solar radiation by the 
negatively charged hydrogen ions is crucial to reach a correct thermal balance at our planet. 

In general, the light absorption by an arbitrary atomic system is determined by separate 
contributions from the corresponding bound-bound, bound-free and free-free transitions in 
this system. As mentioned above the negatively charged hydrogen ion has only one bound 
(ground) 1^5—state, i.e. for the H~ ion there is no contribution from the bound-bound 
transitions. The bound-free transitions in the H“ ion correspond to the photodetachment, 
Eq.(IT]), while the free-free (electron) transitions formally represent the radiation-stimulated 
electron scattering in the held of the neutral hydrogen atom (or inverse bremsstrahlung, for 
short) 

e~ + hv + Y{ = Y{ + e~ (2) 


where the kinetic energy of the hnal electron differs from its incident energy due to the 
photon’s absorption. If the kinetic energy of the both incident photon and electron in Eq. @ 
are small, then the hnal state of the hydrogen atom coincides with its original (ground) state. 
For our Sun the role of inverse bremsstrahlung is substantial for infra-red photons with A ~ 

12,000 - 20,000 A . 

In this study we determine the photodetachment cross-section of the negatively charged 
hydrogen ion (or H” ion). In particular, we want to increse substantially the overall accuracy 
of numerical calculations of the photodetachment cross-section of the H“ ion. Gurrently, the 
photodetachment cross-section of the H“ ion must be determined to very high accuracy (± 
0.1 - 0.2 %). The moderate accuracy known from earlier works (± 2 % - 10 %) is not sufficient 
to describe effects directly related to the climate change and surface temperature variations 
at our planet. In our calculations we apply our highly accurate bound state wave functions 
constructed for the two-electron H~ ion. The crucial part of the method is analytical and 
numerical computations of three-particle integrals which contain spherical Bessel functions 
of the hrst kind. Here we solved this problem and now we can determine a large number of 
different three-particle integrals with Bessel functions jhiPe^) to very high accuracy. This 
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paper opens a new avenne in accnrate nnmerical compntations of the photodetachment 
cross-sections of negatively charged hydrogen ions: (or T“), ^H~ (or D“) and 

By using our current methods we can determine the photodetachment cross-sections 
for all these ions of hydrogen isotopes. Applications of our procedure are not limited to 
the cases of the ground state of the hydrogen atom(s), but also include a number of excited 
states. We also discuss inverse bremsstrahlung, i.e. radiation emitted by an electron which 
moves in the held of the hydrogen atom and a large number of bound state properties of 
the negatively charged hydrogen ions. The computed properies allows us to determine the 
lowest-order QED corrections for these ions. 


II. PHOTODETACHMENT OF THE H ION 


In this Section we derive and discuss the explicit formulas which are used for numerical 
computation of the photodetachment cross-section of the H~ ion. Note that the absorption 
of radiation by the negatively charged hydrogen ions has been investigated in numerous 
earlier studies (see, e.g., ^ and references therein). An extensive review of different 

experimental and computational procedures used to study the interaction of light with the 
hydrogen ions can be found in [^. This paper also contains all essential references 
related to this problem up to 2014. In general, the photodetachment corresponds to the 
bound-free (optical) transitions, since the hnal state of one of the two electrons is an unbound 
state, or, in other words, the state from the unbound spectra of the H” ion. For the two- 
electron ion it is possible to obtain the closed analytical formula for the photodetachment 
cross-section, if the non-relativistic (or dipole) approximation is used. The derivation of this 
formula is considered below in this Section. 

Let us assume that the incident (i.e. second) electron was bound to the neutral hydrogen 
atom with the binding energy e* = —I, where / = Xi is the ionization potential of the H” 
ion. The incident photon has the momentum k and the frequency u (or energy hu = u). In 
the hnal state the emitted photo-electron moves as a free particle with the momentum p and 
energy e/. Since p is a continuous variable, the photodetachment cross-section is written in 
the form 

da = 271 \ Ai^f \'^ 6{-I+ u - (3) 

ATT) 

where Ai^f is the transition amplitude (see, Eq.dH]), below) and the wave function of the 
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final state is normalized per one particle in unit volume V = 1. All wave functions and 
expressions used in computations of the transition amplitude Ai^i are assumed to have a 
proper permutation symmetry upon spin-spatial coordinates of the two electrons 1 and 2. 
This allows one to write the formula, Eq.Q, and all formulas below in one-electron form. 

The delta-function in Eq. ([3]) is excluded by integrating over the momenta p of the photo¬ 
electron. Indeed, by using the formula d^p = p^d | p | do = e | p | dedo and performing 
the integration over de we can remove this delta-function and obtain the following general 
expression for the cross-section 


da = 


e I P 
(2vr)^ 






do 


(4) 


In the non-relativistic approximation {hu <C nieA = Eg ~ e) used below we can replace e by 
rUe- The derivation of the analytical formula for the transition amplitude Ai^f is drastically 
simplihed (see, e.g., 18|]) with the use of transverse gauge for the vector potential A, i.e. 


divA = 0. In this case the transition amplitude Ai^f is written as a scalar product of the 
vector potential A and transition current j*^/, i.e. Ai^f = —eA ■ ji^f- The expressions for 
the vector potential A (in the case when divA = 0) and for the transition current jj^./ can 
be found in standard QED books (see, e.g., [l8]). With these expressions we can write for 
the transition amplitude 


Ai^f — —eA ■ — 



(5) 


where 


Mi^f = J J 2)(e • tti) exp(*k/ ■ ri)V’i(l, 2 )d^rid^r 2 (6) 

and '^i(l,2) and '^/(1,2) are the wave functions of the incident and final atomic systems. 
These two functions are assumed to be properly symmetrized in respect to all possible 
electron-electron permutations and each of these functions has a unit norm. Also, in this 
equation ai are the three Dirac matrixes of the hrst electron, e and kj are two vectors which 
describe polarization and direction of propagation of the incident photon, respectively. 

In this study all electrons are considered as non-relativistic particles, while the energy 
of the incident photon hu{= co) is assumed to be larger than I, but substantially smaller 
than the energy of electron at rest, i.e. hu <C rrieC^. This means that the hnal velocity of 
the photoelectron is small, i.e. we are dealing with the non-relativistic problem. Therefore, 


6 








the Dirac’s a matrixes in Eq.(j6]) can be replaced by the corresponding components of the 
velocity vi = —^Vi. Moreover, in the non-relativistic approximation we can replace the 
factor exp(*kj ■ ri) by uinty. This gives the following formula for the photodetachment 
cross-section 


da,, = 


Pe 


2mpTTUj 




.( 


d 


d 


^ dv2 


Th- 


(ri, r 2 )d^rid^r 2 || do 


(7) 


where ej is the vector which represents the polarization of the incident photon and pe is 
the momentum of the photo-electron. In the case of photodetachment, Eq.(IT]), the notation 
fi{ri, r 2 ) in Eq.(I7]) stands for the wave function of the final state (after photodetachment), 
while TH-(ri,r 2 ) is the wave function of the ground l^S*—state of the hydrogen ion. In 
atomic units the last formula from Eq.([7]) takes the form (compare with 19|, 20| ) 


da,, = aon 


Pe 


271U 


^}*(ri,r2) 


/ d u 
. I -+- 


d 




Th- (ri, r2)d^rid^r2} pdo 


( 8 ) 


^clri dr2' 

2 

where c is the speed of light in vacuum, oq is the Bohr radius and Uq = ^ is the atomic 
velocity. From here one hnds that in atomic units {a.u.) the dimensionless ratio equals 
a~^ ~ 137, where a is the hne structure constant. All other values in Eq.(|8]) must be 
expressed in atomic units. Finally, one hnds for the differential cross-section of the photode- 
tachement 


da,, = —aa\ 

71 




d 


where Pe =\ Pe \ and the factor equals to the ionization potential of the negatively 
charged hydrogen ion H~. In these notations one hnds hu = -1- or u = |pg -f 

in atomic units (it was used in derivation of Eq.(|n])). The formula, Eq.([7]) - Eq.([n]), are used 
in numerical calculations of the diherential and total cross-sections of photodetachment of 
the negatively charged H“ ion. Note that the formula, Eq.([7]), is written in the two forms: 
one-electron form and two-electron form. One-electron form is used in calculations, while 
two-electron form is needed to show the explicit invariance of the hnal expressions upon all 
possible electron-electron permutations. To explain all essential details of such calculations 
we need to derive the expressions which include the approximate wave functions of the bound 
DS'—state in the H“ ion(s) and wave functions of the hnal states, i.e. states which arise 
after photodetachment. Diherent forms of the hnal state wave functions explain diherent 
formulas which are used in the modern literature for these cross-sections. 
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III. WAVE FUNCTIONS OF THE NEGATIVELY CHARGED HYDROGEN IONS 


The wave functions of the ground states in the hydrogen ions (or T“), 

^H~ (or D“) and ^H“) are approximated with the use of different variational expansions. One 
of the best such expansions is the exponential variational expansion in relative/perimetric 
coordinates coordinates. For the ground S'—state of the H~ ion it takes the form 2l| 


1 


N 


TH-(ri,r2) = exp(-air32 - Argi - 7ir2i)+exp(-/3ir32 - Oirgi-7ir2i) (10) 

V ^ i=l 

where A 2 is the permutation of the two identical particles 1 and 2 (electrons), the subscript 
3 designates the hydrogen nucleus. In this notation r 2 i =| r 2 —ri | is the electron-electron (or 
correlation) coordinate and r^i =\ r 3 —r* | are the corresponding electron-nuclear coordinates 
{i = 1,2). Each of these three relative coordinates ^ 32 , r 3 i, r 2 i is a scalar which is rotationally 
and translationally invariant. Also in this formula Ci {i = 1, ..., N) are the linear variational 
coefficients of the approximate (or trial) wave function. The 3V—parameters (i = 

1 , ..., N) are the non-linear parameters (or varied parameters) of the trial wave function, 
Eq. lfTOi) . As follows from the results of our earlier studies (see, e.g., 0 . H) the exponential 
variational expansion, Eq. flTnil . provides very high numerical accuracy for any bound state 
in arbitrary Coulomb three-body systems, i.e. for systems with arbitrary particle masses 
and electric charges. For the H“ ions such an accuracy is much higher than accuracy which 
any other variational expansion may provide, if the same number of basis functions N is 
used. Note also that from Eq. ffTOD one finds the following formula 


f) F) ^ ^ 

e/ ■ (j— + = (e/ • + U, ■ £i) 


( 11 ) 


h' ■ “■) [flr:' 


r ^ . 

( \ 

r ^ . 


+ (,6/ ■n2j 

0 — 

2 


’■ 1 ' a 

= le/ ■ ni ' 


where n,- = ^ is the unit vector of the Ath electron and z = 1, 2. Note that in this 

expression all terms which contain the 7 ^ parameters cancel each other, since we always 
dr-2.1 


have TpJ' = 

ari 


dV2 


The presence of a large number of the varied non-linear parameters aj,/Sj, 7 i in Eq. flTOl) 
allows one to construct very compact, extremely flexible and highly accurate approxima¬ 
tions to the actual wave functions of arbitrary three-body systems, including the negatively 
charged hydrogen ions. The results of our calculations obtained with the use of the exponen¬ 
tial expansion Eq. flTOl) . can be found in Table I. This Table contains only the total energies 
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(in atomic units) and ionization potentials xi (iii eV). In addition to the values shown 
in Table I we have determined a large number of expectation values which can be used to 
evaluate various bound state properties and probabilities of different processes/reactions in 
these ions. Some of these expectation values determined with our wave functions are shown 
in Tables II and III. The meaning of the notations used to designate some bound state prop¬ 
erties is explained in j^|. Table II illustrate actual convergence rates for different properties 
upon the total number of basis functions used. It allows us to predict the number of stable 
decimal digits in these and other similar bound state properties. Table III contains a large 
number of expectation values (only stable decimal digits are shown) determined in highly 
accurate computations. Note that the total non-relativistic energies and other expectation 
values obtained in our calculations are the most accurate values ever obtained for isotope 
substituted hydrogen ions: (or T“), (or D“) and ^H“. Convergence of these expec¬ 

tation values can be understood from Table IV. These expectation values indicate the overall 
quality of our wave functions used in calculations of the photodetachment cross-section of 
the H“ ion(s). Another well known test for bound state wave functions of Coulomb few- 
□ody systems is numerical coincidence of the predicted and computed cusp values (see, e.g., 


2l|). Recently, we have found that the quality of bound state wave functions can be tested 


in accurate evaluations of the lowest-order QED correction for the hydrogen ions (see Ap¬ 
pendix I). Our values for the lowest-order QED corrections determined to high accuracy for 
the °°H“, (or T“), (or D“) and ions are presentd in Table V. 

It should be mentioned that in our calculations of the and hydrogen ions 

we have used the following updated values for the nuclear masses of the tritium, deuterium 
and protium (in MeV/(?) 


me = 0.510998910 , = 938.272046 

md = 1875.612859 , mt = 2808.920906 


( 12 ) 


In our calculations these masses are considered as exact. The masses of these three hydro- 
genic nuclei have recently been determined in high-energy experiments to better accuracy 
than they were known in the middle of 1990’s. Usually, these masses are expressed in special 
high-energy mass units {MeV/c^). Numerical value of the electron mass mg (or mgC^) can 
be used to re-calculate these masses to atomic units. Highly accurate computations of the 
ground states in the (or T“), (or D“) and ^H~ ions with these nuclear masses have 
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never been performed. 


IV. FINAL STATE WAVE FUNCTIONS 


After photodetachment of the H“ ion, Eq.(IT]), in the hnal state we have the neutral 
hydrogen atom in one of its bound states and ‘free’ electron which moves in the field of 
the hydrogen atom with the kinetic energy E^. = where h = ^ is the reduced Planck 
constant (or Dirac constant), k is the wave number and rrie is the electron mass. The 
corresponding wave function of the final state is represented as a product of the bound state 
wave function of the hydrogen atom and the wave function of the free electron. The wave 
unction of the bound (n£m)—state of the hydrogen atom is written in the form (see, e.g.. 


22[) 0,0) = i?„^(Q,r)Y£,n(0,0), where Y£m(0,0) = >£m(n) is a spherical harmonic 

and Rne{Q,r) is the one-electron radial function. The radial function Rne{Q,r) takes the 


form (see, e.g., [22|) 


Rnii^Q ) ^) 


Q{n-^- 1)! r 2Qr iUi f (-1) 


rn\ {nRtj\ 


L n J 


n—i—l 

i E 

A:=0 


k\ 



\2Qr 


L n 


exp 


]x 

(-?) (-) 


where Q (= 1) is the nuclear charge, while n and £ are the quantum numbers of this bound 
state. Note that all radial functions defined by Eq. fllSp have unit norms. The six following 
functions of the hydrogen atom are very important for calculations of the photodetachment 
cross-sections: 

() . A2i(r)= ^ 


Rio{r) = 2exp(-r) 
Rzo{r) = 


R2o{r) = ^exp(-^)(l 


2 / 2r 2r‘^\ . r. 

—;=(!-1-) exp( — 

3\/3'^ 3 27^ 3' 


2' 

Rsiir) = 


^rexp(—) (14) 

2x/6 2^ ^ ^ 


8 Yl _ !() exp(-!:) 


27V6^~ 6' 

The final states with these hydrogen wave functions provide largest numerical contributions 
in the photodetachment cross-section. Note also that in this study all angular parts of these 
hydrogenic functions (spherical harmonics) are included in the formulas for matrix elements 
(see, e.g., Eq. llTSll below). 

The wave function of the free electron is represented as the sum of the plane wave 0(r) = 
exp(?k ■ r) and converging spherical waves which are usually designated as 0|,“^(r) 


22|. In 


the non-relativistic dipole approximation used here the following selection rule is applied: 
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transitions from the incident S{L = 0)—state of the ion can only be into the hnal 
P{L = 1)—state of the two-body (hnal) system: the H-atom pins free electron. It follows 
from here that in the hnal states we can detect the H-atom in one of its s{i = 0)—states 
and free electron which moves away in the p-wave. This is the sp—channel of the 
ion photodetachment which prodnces the largest contribntion into photodetachment cross- 
section of the ion. Another possibility is the formation of the hnal H-atom in one of its 
three p{i = 1)—states and free electron which moves away in the s—wave. This represents 
the ps—channel in the photodetachment of the H“ ion. 


The wave fnnction of the hnal state for the sp—channel takes the form (see, e.g., 18|) 


^ RnoiQ ^ l,J’3l)(ne ■ n2)<I>pi(r32) ^ (lie ■ ll2)i?„o(l’3l)jl(Pel’32) 


(15) 


where no = — is the nnit vector of the second (nnbonnd) electron and m is the 

direction of the momentnm of the ontgoing photo-electron pe, i.e. the second electron 
which becomes nnbonnd after the photodetachment of the hydrogen negatively charged 
ion. In Section HI this vector has been designated by p, i.e. pe = p. The notation 
Rne=o{Q = = -Rno(’" 3 i) stands for the radial fnnction of the bonnd state of the 

hydrogen atom with qnantnm nnmbers n > 1 and i = 0. In this eqnation the notation 
*^^ 61 (^ 32 ) = Ji(Pe^ 32 ) stands for the wave fnnction of the electron which becomes ‘free’ after 
photodetachment. Here and below ji{x) and jo{x) are the spherical Bessel fnnctions of the 
hrst and zero-orders, respectively. Snch a wave fnnction represents an electron which moves 
in the held of the nentral hydrogen atom. This electron has an angnlar momentnm ^ = 1 
and absolnte momentnm pe =| Pe |. 

The wave fnnction of the hnal state for the ps—channel can be written in the following 
form 


^/>/(l, 2) = Rni=i{Q = 1, r3i)Yim(ni)^<I)p£=o(?^32) = • Rni{r‘ii)io{per‘i2) (16) 

^Pe ^Pe 

where Rni{Q = l,'r 3 i) = Rnii^si) is the radial fnnction of the bonnd | nl)-state of 
the hydrogen atom with the principal qnantnm nnmbern {n > 2) and angnlar qnan¬ 
tnm nnmber £ = 1 . Also, in Eq. flTB]) the notation Yim(x) stands for the correspond¬ 
ing spherical harmonics 2^, while the vector ni = ^ has the following components: 
(ni)+i = -^{nx + tUy), (ni)o = n^, (ni)+i = -^{ux - tUy), where nx,ny,n^ are three 
Cartesian components of the nnit vector rii. Note that all calculations for the ps—channel 
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in the photodetachment of the H~ ion are signihcantly more complicated than for the 
sp—channel. One obvious complication can directly be seen from Eq. flT6|) . where the £- 
nal wave function is a vector. In this study we restrict ourselves to the consideraion of 
the ‘traditional’ sp—channel only. It is clear a priori that the overall contributions of all 
ps—channels into photodetachment of the H~ ions in actual stars is very small. Excitation 
energies which are needed to form the hnal hydrogen atoms in the np-states are extremely 
high (AE > 10.2043 eV). The corresponding ‘equilibrium’ temperatures at the stellar sur¬ 
face are above 10^ K (they correspond to very hot 0-stars also known as the Wolf-Rayet 
stars). However, as mentioned in the Introduction for such stars photodetachment of the 
H“ ions dos not play any role. 


V. FORMULAS FOR THE MATRIX ELEMENTS 


In this Section we derive analytical expressions for the photodetachment cross-section(s). 
First, let us discuss the angular integration in the expression for Mj^/ included into Eq.(IH]). 
Below we restrict ourselves to the case of the sp—channel only. In this case the hnal state 
wave function is represented in the form of Eq. flT5|l and the formula for the amplitude 

is written in the following form 


- 


A(pe • n32)(e/ ■ n3i)Ri -h A(pe • n32)(e/ ■ n32)R2 


dr3i dr.q 


31«i32 
2 


(17) 


■ n32)(e/ • n3i)Bi + A(p^ ■ ■ n32)B2^r.^idr3ir.^2‘l'‘^32dn„dn,2 

where A, Bi and B 2 are some scalar expressions (see below). Result of angular integration of 
the hrst term in the last equation equals zero identically, since its integrand is a scalar upon 
n 3 i. The integration of the second term produces the factor ^(e/ ■ Ue). By introducing the 
explicit expressions for the A, Bi and B 2 scalar values we obtain the following formula for 
the transition amplitude Mi^f\ 


N 

Mi^f = 47r?(ne ■ ej) ^ Cj 
i=i 


f +00 /■+00 /■r3i+r32 


ad 


' k31-7’32| 


Rno{r3l)jl{Per32) X 


(18) 


exp{-ajr32 - PjTsi - 'yjr2i)r32r3ir2idr32dr3idr2i + fdj x 

^+00 ^+00 j■rzl+rz2 

/ / / Rno(?^3i)ji(Per32)exp(-/3jr32 - ajr3i-7jr2i)r32r3ir2idr32dr3idr2i 

JO JO J|)'3i-r32| 

where the notation ‘h£=i(pe;r 2 ) designates the wave function of the ‘free’ electron which 
moves in the held of the neutral hydrogen atom. As one can see from Eq. ffTSj) numerical 
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calculations of the transition amplitude Mi^j are reduced to the calculation of the three- 
body integrals in relative coordinates r 32 , rai, r 2 i. Calculations of such integrals are discussed 
in detail in the next Section. Here we want to note that the photodetachment cross-section 
for the sp—channel always contains the factor | p which is proportional to the factor 

(Pe • e)^ = Pg(ne • e)^ (here and below e = ej). 

To determine the factor (pe • e)^ = p^Ue ■ e)^ we note that in 99.99 % of all stellar 
photospheres the hydrogen ions are free-oriented. This means that the expression from 
Eq. flT7|) must be averaged over possible directions of the unit vectors ei and 62 which describe 
the polarization of the incident photon. In the incident light wave we always have ei T e 2 T 
kf, where the unit vector kj corresponds to the direction of propagation of this light wave, or 
the corresponding light quantum (photon). An arbitrary scalar expression can be averaged 
over all possible directions of the unit vectors ei and 62 with the use of the following formulas 


2 f')f') k') 


(19) 


and 

(a ■ e)(b ■ e) = ^[(a ■ b) - (a ■ kf){h ■ kj)] = i(a x kj)(b x kf) (20) 

where a and b are the two arbitrary vectors and the notation x x y stands for the vector 
product of the two vectors x and y. From Eq. (|2nll one finds that 

(ng-e)2 = i(ne X kj)2 (21) 

Now, the explicit formula for the differential cross-section, Eq.(jn]), takes the form 


da^ 


SiraalJCi )(ne 

^Pe+l ^ 


X k/)' 






do 


( 22 ) 


where K, is some additonal (numerical) factor and the matrix element (photodetach¬ 

ment amplitude) is written in the form 




N 

i=i 


r+oo I'-l-oo I•r3l+r32 


C^j I I I Rno{r3l)jl{Per32) X 

JO JO 2|r3i-r32| 

exp(-Q;jr32 - PjVsi - 7i’^2i)’^32’^3H2idr32dr3i(ir2i + (3j x (23) 

^+00 r+oo rr 3 i+r 32 

/ / / Rno{r3i)ji{Per32)exp{-(]jr32-ajr3i-'Jjr2i)r32r3ir2idr32dr3idr2i 

Jo Jo -tlrsi—r 32 | 


The factor /C is responsible for missing normalization constants in the incident and final wave 
functions. The same factor includes one part of the hnite mass corrections for the protium. 
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deuterium and tritium (hydrogen isotopes). Another part of such a correction comes from 
the photodetachment amplitude 

As follows from Eq. fl22p the photodetachment cross-section is always proportional to the 
factor (ne x = sin^ 0, where 0 is the angle between the directions of propagation of 
the incident photon kj and hnal electron He (or photo-electron). This results corresponds 
to the well known dipole approximation. As follows from Eq. fl2^ the matrix element 
does not depend upon any of the angular variables. Therefore, the formula for the total 
cross-section of the photodetachment takes the form 

647r^ 


(Tr, = 


-aa'i 


^Pi + T^ 






(24) 


Now, we need to perform the last step of our procedure and calculate three-electron integrals 
in the relative coordinates. Note that each of these integrals contains one spherical Bessel 
function jq{pr 32 ), be. for the sp—channel and jo{prz 2 ) for the ps—channel of the H~ 

photodetachment. 


VI. THREE-BODY INTEGRALS IN RELATIVE COORDINATES 


As we have shown above the problem of analytical and numerical computations of matrix 
elements arising in Eq.([7]) and Eq. lfTSj) with the trial wave function dehned in Eq. flT0|) is 
reduced to the computation of the following three-body integrals with the Bessel functions 
Jo(^^2) = jo{p^ 2 ) for the ps—channel and ji{kr 2 ) = Ji(pr’2) for the sp—channel 241 

^+00 r+00 rr 3 i+r 32 

/ / / 0n£(r3i)jg(pr32)exp(-a;r32 -/^rsi - 7r2i)r32r3ir2idr32dr3i(Jr2i (25) 

Jo Jo J|)'3i-r32| 

where g = 0, 1 and r2i =| r2 — ri | is the electron-electron coordinate which varies between 
the following limits | r32 — r3i |< r2i < r32 -|- r^i. In Eq. fl25]) the notation stands 

for the unit norm wave functions i?„o(^3i) and Rniir^i) of the bound states of the hydrogen 
atom, Eqs. llT^ and (ITTll . All other notations in Eq. (|25|l are dehned above in Eq. llTOil . The 
integrals dehned by Eq. fl25|l are computed analytically in perimetric coordinates Ui,U 2 ,U 3 , 
where Ui = |(r3i -f rsi - r32),M2 = |(?^32 + ?^2i - r^i) and M3 = |(r3i -F rs 2 - ^21). The 
three perimetric coordinates ui,U 2 , M3 are independent of each other and each of them varies 
between 0 and -|-cx). 

Let us derive the closed analytical expression for the integrals, Eq. fl25]l with g = 0 and 


14 







g = 1. First, we need to obtain the general formula for the closely related auxiliary three- 
body integral /9, 7 ) which is written in the form 

^+00 ^+00 ^r3i+r32 

rfc,i,n(a,/?,7) = / / / ’’ 32 ^ 3 i’^ 2 i exp(-ar32 -/^rsi - 7r2i)cir32dr3i(ir2i (26) 

Jo Jo 2|r3i—r32| 

where fc > 0 , / > 0 , n > 0 and a-l-/ 5 > 0 ,a-|- 7>0 and /? -f 7 > 0 (see below). Analytical 
computation of this integral has extensively been explained in a number of erlier studies. 
Correspondingly, below we restrict ourselves only to a few following remarks. In perimetric 
coordinates the integral, Eq. fj26ll . takes the form 


Tk,l,n{(y,P,7) 



(a + / 3 )m 3 - (a + 7 )m 2 - (/? + 7 )«i] (“2 + 

(ui + U'i)\ui + U 2 Y duidu2du‘i (27) 


where we took into account the fact that the Jacobian of transformation from the relative 
(^" 32 , ^ 31 .^ 21 ) to perimetric coordinates (mi,M 2 ,M 3 ) equals 2. The integration over three inde¬ 
pendent perimetric coordinates u, (0 < Uj < cxo) in Eq. (l271) is simple and explicit formula 
for the ;(a,/3, 7 ) integral is reduced to the form 

{l-li + ki)\ {k-ki + ni)\ {n-ni + li)\ 


U;,;„(a,/3,7) = 2 E E E Ct'Cl'Cl'- 


ki =0 li =0 ni =0 
I n 


Of + 


/ — 


a+ 7 ) 


/c—/ci+ni + l 


{(d + l) 


n—ni+^i +1 


ni+fci /i+ni 


= 2 • A’’ • /’ • 77^ V V V _ 

I^Oh^Or^O (“ + /5)^“'l+''l+^(a + -y)fc-fcl+"-l + l(/5 + ^'jn-m+h+l 


(28) 


where C™ are the binomial coefficients (= number of combinations from k by m) (see, e.g.. 


25|h In many cases all three integer numbers k,l,m are relatively small (< 3) and this 


simplihes numerical applications of the formula, Eq.( 

The parameters a,/3 and 7 in Eqs. fl2^ - (|28|) can be arbitrary real and/or complex 
numbers. Furthermore, one of these three parameters can be equal zero, but the three 
principal conditions for these parameters a-|-/J> 0 ,a-l- 7>0 and /3 -|- 7 > 0 must always 
be obeyed. If these three parameters are complex, then these three conditions are written 
in the form Re{a - 1 - /?) > 0, Re{a - 1 - 7 ) > 0 and Re{(3 - 1 - 7 ) > 0. These conditions are needed 
to guarantee convergence of the auxiliary three-body integrals Ffc;;;„(a,/J, 7 ), Eq. fl2^ . 

Now, we can derive the explicit formulas for the analogous three-body integrals which 
contain the Bessel functions jo(p^ 2 ) and/or ji(pr 2 ), where p = Pe- First, let us obtain the 
computational formula for the integral which contains the jo{p^o 2 ) function 


’ ’ Jo 


+00 /■-l-oo I•r32+r3l 


'732—r3l| 


’^32’^31^2iJo(P?^32) X 


exp(-ar32 - (droi - 7’^2i)dr32dr3idr2i 


(29) 
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By using the formula jo{x) = one finds 

.l)Qp2g 

To. 

q=0 






-i)V 


<lrTl 


Qmax 

~ § (29+1)! 


rfc+2g;Z;n(a,/5,7) (30) 


where rfc.;.„(a,/3, 7 ) is the integral defined in Eq. fimi . 

The integral l3, 'y; p) in the last equation converges for relatively small p, but for 

p < 1 it converges very rapidly. In reality, the maximal values of the index q (i.e. qmax) 
in Eq. fl30|) are always hnite. Numerical investigations indicate that to stabilize 15 decimal 
digits for fc < 1 one needs to use in Eq. fl30|) qmax = 5 and qmax = 10. For p > 2 the value 
of qmax rapidly increases up to 40 - 50. The same conclusion is true about the convergence 
of the three-body integrals with the Bessel function ji{x) = This integral takes 

the form 

, . r +00 I'+oo |•r32+r3l 

= / / / r32r3ir2iji(pr32)exp(-ar32-/^rsi-yraOdrsadrgidrai 

Jo Jo J\r32-r3i\ 

-iy{2q + 2)p^^+\ 


10 

00 

= E' 

q=0 


"T }iJ^2q+l\l-,niS^ 1 y I 7 ) 


(31) 


(2g + 3)! 

The knowledge of the three-body integrals with the Bessel functions io{x) and ii{x) al¬ 
low one to determine all integrals, Eq. fl2^ . needed in calculations of the photodetachment 
cross-section of the negatively charged hydrogen ion. The integrals j3^'y]p) are 

needed to calculate the photodetachment cross-sections in the sp—channel, while the inte¬ 
grals /3,'y,p) are needed to evaluate analogous cross-sections in the ps—channel. 

Finally, let us present the formulas for the three-body integrals with the spher¬ 
ical Bessel functions jo{Pe‘>^ 3 i) and ji(per’ 32 ) and a few radial hydrogenic functions 
Rio{r), R 2 o{r), R 2 i{r), Rsiir), where 


Rio{r) = 2Q^/Qexp{-Qr) , i? 2 o(?^) = - ^Qr)exp(-^) 




R2iir) = ^^Qrexp{-^) 


^/2 

Rsi{r) = (l - ^Qr)ea:p(-^) 


27\/6 6 ' 

where Q = 1 for the hyrogen atom. For the photodetachment of the hydrogen negative 
ion only the integrals which contain the following products i?io(’" 3 i)ji(p^" 32 ), -R 2 o(’" 3 i)ji(p^ 32 ) 


and 7?2i(?"3i)jo(p^"32), .R 3 i(?" 3 i)jo(p^" 32 ) are important. As follows from Eqs. 
integrals are 

2Q^jQB[]lya,(3 + Q,-f;p) , 


QVQ ^(1) ( Q , Q ^ 

5 l;l;l («,/5 + Y, 7 ;P) 


Q\/Qd(i) ( q , q ^ 

-^—B\.yya,y + —,7;p) 


and these 

(32) 

(33) 
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and 


Q^VQ d(o) o , Q 
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5};2;i(a,/? + —,7;p) 


8Q\/Q (0) / Q , Q ^ 4Q^i/Q (0) ^ ^ , Q ^ 

(a. + -. Ti p) - ^^Bp,p (a. 7; P) 


(34) 

(35) 


respectively. Analogous formulas can be written for the integrals which include other radial 
hydrogenic functions, e.g., Rno{r) and Rmiir)^ where n > 3 and m > 4. However, the 
overall contribution of such integrals to the photodetachment cross-section of the H“ ion 
(Q = 1 in all formulas above) is very small and they can be neglected in the lowest order 
approximation. 

To test our formulas we perform calculations of the photodetachment cross-sections of 
the °°H“, (or T“), (or D“) and ions. These cross-sectons are determined to 
high accuracy and for a large number of photo-electron moment a/energies (see Tables VI - 
VIII). It is shown in 2^ that the additional factor /C in Eqs. fl22l) and fl241) equals ^ which 
is the normalization constant of the angular part of the three-particle (or two-electron) wave 
functions of the H“ ions. This leads to the following formula for the total cross-section of 
photodetachment 


_ 8 2 / Pe \ 

— Q (,2 1 2 j 






(36) 


for the °°H ion and 




:aal(l + 


. 2 r 


Pe 


Ri 




(37) 


3--''V- ■ M+n Lp2(l+,^) + ^,2. 

for the hydrogen isotopes with finite nuclear masses. In Eq. fl37p the notation M stands for the 
mass of the nucleus of the corresponding hydrogen isotope (in atomic units). Results of our 
computations of the photodetachment cross-sections of the hydrogen negatively charged ions 
can be found in Tables VI, VII and VIII. In these calculations we have used our variational 
wave functions with N = 300 and N = 350 for the ion and N = 350 for other 

hydrogen isotopes {N is the total number of basis functions used). The corresponding 
total energies for each trial wave function (in atomic units) are given in Table VI. For the 
ion our photodetachment cross-section is comparable with the analugous cross-section 


from 


16| (only our cross-section is more accurate). Comparison with other papers where 


the photodetachment cross-section of the negatively charged hydrogen °°H ion has been 
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determined in the velocity representation (see references in [1^) indicate clearly that our 
photodetachment cross-section (for this ion) are in good agreement with the cross-sections 
known from earlier works. Note again that the °°H~ ion is a model atomic (or two-electron) 
system with the inhnite nuclear mass. To the best of our knowledge the photodetachment 
cross-sections of the hydrogen ions with hnite nuclear masses have never been determined 
in earlier works. In reality photodetachment of the negatively charged hydrogen ions with 
the hnite nuclear masses requires an additional consideration. In our current approach we 
just introduced the corresponding mass-dependent corrections in the normalization factors 
of the hydrogenic wave functions (for the hnal hydrogen atom). All numerical calculations 
have been performed by using the formulas Eqs. fl36p - fl37p . A large number of additional 
results for the photodetachment cross-sections of the negatively charged hydrogen ions will 

29l |. which inlcudes the photodetachment cross-sections in 


be presented in our next paper 


those cases when the hnal hydrogen atom is formed in the excited s—states. 


VII. A SIMPLE METHOD TO EVALUATE CROSS-SECTIONI 


The photodetachment cross-section of the negatively charged II~ ion in the sp—channel 
can approximately be evaluated with the use of the following asymptotic method. Let us 
designate the total number of bound electrons in atomic system by Np, while the notation Q 
stands below for the nuclear charge Q. As is well known (see, e.g., 13| - 0 and references 
therein) the long range, single-electron asymptotic of the wave function of A^^e—electron 
atomic system with the nuclear charge Q is written in the form 


I T(r) \=Crt ^exp(—tr) (38) 

where C is some numerical constant, Z = Q — N^ + l, while t = ^/2J^ and A = Xi is the hrst 
ionization potential which corresponds to the dissociation process H” = H -|- e“. For the 
ion one hnds in Eq. (l3^ : Q = 1, Ng = 2, and therefore, Z = 0 and the long range asymptotic 
of the wave function, Eq. fl38l) . is represented in the Yukawa-type form | \k(r) |= ^exp(—tr), 
where C is a constant which must provide the best correspondence with the highly accurate 
wave function of the II~ ion at large r. The highly accurate wave function of the II~ ion is 
assumed to be known. The explicit formula which is used to determine the constant C is 
C =1 \['H-(r) I rexp(tr), where \I'H-(r) = TH-(r 32 = r,rsi = 0,r2i = r). As follows from 
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the results of actual computations the factor C in Eq. fl38|) does not change for relatively 
large interval of variations of r, e.g., for r bounded between 7.5 and 35 atomic units. The 
determined numerical value of this constant C for the negatively charged hydrogen ions 
varies between 0.3562404 (°°H~) and ~ 0.3565232 (^H~). These values are larger than 
the corresponding values obtained in [l7| by a factor 2\/2. It is important to note that 
this simple asymptotic method is essentially a single-electron approach. Therefore, only 
three Cartesian coordinates r and one radial variable r are used below. These coordinates 
correspond to the electron which was originally bound to the neutral hydrogen atom. At 
the hnal stage of the process this electron becomes free. 

The photodetachment cross-section a{H~) of the H“ ion is written in the form which is 
similar to Eq.(|9]) (see, e.g., jl2| ) 


da(R^) = I ^ ““°2^ ^ 

The normalized wave function of the incident state is exp(—tr) = ^exp(—tr), while 

the wave function of the hnal state is exp(*pe • r). By substituting these exressions into 
Eq. fl3^ one hnds 


d 


/ 1 \ 1 /*H"CXD 

da(H~) = aal(— -^)(e/ ■ Uef- \ — / 'ijji{r) sm{per)rdr do 

^ JiTTUJ ' J)(s Jo 

■ Svrpf 


Pe Jo 
+ 00 


2 oTrPg , ( \ ' ( \j| 2 j 

aa^ - ief • ng) • — / exp[—'yr)sm[per)dr do 

CO Vp Jo 


(40) 


where Pe = Pei^e and Ue is the unit vector which determines the direction of propagation of 
the emitted photo-electron. Also, in this equation here a is the hne structure constant and 
Oo is the Bohr radius and 7 = t = 1 / 2 / 1 . A few steps of additional transformations lead to 
the following hnal formula 


(icr(H ) = (k/- x nei^do (41) 

^Pe + 7 ^ 

where we used the identity (e/ ■ Ue)^ = |(kj x rie)^ and the fact that the energy conservation 
law (for the photon) is written in the form: hu = or in atomic units: u = ^ + \- 

To transform the formula from Eg. (1391) to its hnal form we have chosen the incident (one- 
electron!) wave function in the form of Eq. (l38|) . while the hnal wave function "0/ is written in 
the form ea;p(*Pe ■ r). Also in this equation pe =| Pe I is the momentum of the photo-electron 
(in atomic units) and 7 = 0.23558869473885 (see Table I). 
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By integrating over angular variables {do = sin 6d6d(j)) one finds the total cross-section 
of the photodetachment 


Rzl-Tr^ 

ct(H“) = —:—aal ■ 


Pe 


= 33.9954673 ■ 10"^® • C 


Pe 


cm 


(42) 


3 “ (pl + 'y^)^ ■ (Pl + 'l^)^ 

where C is the constant from Eq. fl38p . This formula with the known value of C can directly 
be used in calculations of the photodetachment cross-section a{H~). This formula essentially 
coincides with Eq.(8) from ll(]. In our case, C = 0.3562404, and therefore, from the last 
formula one finds 


ct(H-) = 4.31427025 ■ 10”^® 


Pe 


{pI + 7 ^)' 


cm^ 


(43) 


This simple method predicts the maximum of the photodetachment cross-section arpn^r. ~ 
41.24335828-10“^® cm^, which is close to the theoretically ‘exact’ maximal value 1^, but 


slightly above it. 

In actual calculations the overall accuracy of the simple formula, Eq. (l4^ . is outstanding. 
This remarkable fact follows, in part, from the comparison of the sizes of atomic system a 
with the wavelengths A of the photons which produce photodetachment of the H“ ion(s). 
It is clear that a ~ oq ~ 0.529A where oq is the Bohr radius and a is the effective radius 
of the atomic H“ ion. On the other hand, the wavelengths A of the photons which produce 
the photodetachment is A ~ 1 • 10®A. This means that A S> a, or, in other words, the 
ratio ^ ~ 2 ■ 10® is a very large number. Briefly, we can say that the incident photon sees 
only an outer-most asymptotics of the three-body wave function. All details of the wave 
function of the H“ ion at small electron-nucleus distances are not important to describe the 
photodetachment of this ion. 

To conclude this Section we have to note that some important details of the photodetach¬ 
ment of the H“ ion cannot be investigated by applying this simple method only. For instance, 
it cannot be used, in principle, to evaluate the transition probabilities for the ps—channel(s) 
and/or predict the final probabilities of the formation of hydrogen atoms in the excited 
s—states for the sp—channels. Nevertheless, this method is often used in astrophysics for 
fast evaluation of the absorbtion coefficient of the H“ ion at different conditions. 
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VIII. INVERSE BREMSSTRAHLUNG IN THE FIELD OF HYDROGEN ATOM 


The knowledge of the photodetachment cross-sections for the sp—channels allows one 
to evaluate the absorbtion of infrared and visible radiation by the H~ ion in photospheres 
of many stars with surface temperatures Tg < 8,250 K. Additional contribution from the 
ps—channels increases the overall accuracy of such evaluations. However, at large wave¬ 
lengths (A > 12, OOOA) a different mechanism begin to play a signihcant role in the absorb¬ 
tion of infrared radiation. Radiation at these wavelengths can be absorbed by ‘free’ electrons 
which move in the fileds of hydrogen atoms. In reality, such an absorbtion by ‘free’ electrons 
lead to the disappearance of a noticeable amount of propagating infrared radiation with 
A > 12, OOOA. It cannot be explained by the photodetachment of the H“ ions. Therefore, 
in this Section we need to discuss the electron scattering in the field of the heavy hydrogen 
atom, Eq.(|2]) and radiation transfer which arise during transitions between the ‘unbound’ 
states of moving electrons. Note that during such a process the energy of solar radiation 
is transferred, in part, to the scattered electrons. Electrons are accelerated, but the mean 
energy of the radiation held decreases. In classical physics this process is called inverse 
bremsstrahlung. Inverse bremsstrahlung is important in many areas of modern physics and 
various applications. In stellar photospheres of hydrogen stars inverse bremsstrahlung is 
resp onsible for absorption of infrared radiation with A > 12,000 A. This includes our Sun 

H- 

As mentioned above at small energies of the incident photon hu < 5eV the original 
(ground) state of the hydrogen atom does not change. This means that for photons of small 
energies we can ignore all possible excitations of the hydrogen atom, i.e. atomic transitions 
from the ground state of the H atom into other excited states. The incident photon is 
absorbed exclusively by the electron which moves in the field of the neutral hydrogen atom. 
In this case the electron scattering at the hydrogen atom (H -|- e“ -|- hz/ = H -|- e“(*)) can 
be considered as a regular (but non-elastic!) scattering of an electron in the held of the 
heavy hydrogen atom H which is electrically neutral. In the lowest order approximation we 
can replace the actual hydrogen atom by its polarization potential V{r) = —where a is 
a positive constant. Note that the potential V{r) is a central and local potential. For the 
°°H atom the numerical value of a in atomic units equals a = |. This means that we need 
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to consider the electron scattering in the held of the following potential 


V(r; a, b) = 


( 44 ) 


(r + by 

which analytically depends upon the two real parameters a and b. It is clear that we are 
dealing with the case of ‘polarization’ inverse bremsstrahlun g, w hich is substantially different 
from the regular (or Coulomb) bremsstrahlung (see, e.g., 27|]). Our goal is to determine 
the cross-section of inverse bremsstrahlung for an electron which moves in the held of the 
potential V (r; a, b), Eg. 044 p . By using diherent numerical values for the a and b parameters 
we can obtain a variety of such cross-sections. The main interest, however, is related to the 
cases when in Eg. 044|) a ~ 2.25 and 6 ~ 0, respectively. 

To calculate the cross-section of inverse bremsstrahlung in the held of the neutral hydro¬ 
gen atom we can use a number of diherent methods which can be separated into two large 
groups. The hrst group of methods is based on the use of very approximate electronic wave 
functions which do not describe electron-electron correlations. However, by using these wave 
functions we can solve ‘exactly’ the corresponding scattering problem. Moreover, by varying 
a number of non-linear parameters included in the incident and hnal wave functions one can 
obtain a good agreement with the known experimental data. The second group of methods 
is based on an explicit construction of the two-electron wave functions. In this approach all 
possible electron correlations are included into wave functions. For bound states such wave 
functions can be determined to very high numerical accuracy. However, it is very difficult to 
obtain accurate numerical solutions of the general scattering problem with truly correlated 
wave functions. The source of the problems is well known, since, e.g., the total number of 
independent variables is different for the incident and hnal wave functions. In this study 
we develop an improved version of single-electron approach. In our approach the hydrogen 
atom is replaced by an ‘effective’ central potential. This interaction potential explicitly 
depends upon a few parameters which can be varied to achieve a better approximation to 
actual scattering data. 

Let us describe the method which can be used to calculate cross-section of inverse 
bremsstrahlung in the held of the hydrogen atom. First, we need to construct the inci¬ 
dent and hnal wave functions which are the stationary one-electron wave functions of the 
continuous spectrum of the radial Schrodinger eguation (see, e.g., |22|) written for one elec¬ 
tron which moves in the held of the central potential, Eg. fl44|) . The unknown wave function 


22 





^'(r) can be represented in the form '^g{r) = Ai{r) exp(z(5£(r)), where Ag{r) is the amplitnde, 
while S£{r) is the phase of the wave fnnction. Note that in the held of any central potential 
the angular momentum £ is a conserving quantum number (or good quantum number). It 
is clear that it should be an additional relation/equation between the amplitude A^lr) and 
a phase Silr) of the wave function, since we cannot simply replace one unknown function by 
two new functions. By using this ‘additional’ equation one reduces the original problem to 
the solution of the following differential equation 


dSjjr) 

dr 


2V (r; a, b) 
k 


it{kr) cos Si{r) — ni{kr) sin ^^(r) 


(45) 


where the phase ^^(r) of the wave function is assumed to be explicitly r—dependent. The 
initial condition for the phase is written in the form 6i{r = 0). With this initial condition 
Eq. (H5ll is reduced to the form of the following integral equation 


de{r) = -j:J^ y{s-,a,b) je{ks) cos6£{s) - ni{ks) sin 6i{s) 


n2 


ds 


(46) 


The equations Eq. (l45ll and/or Eq. (l46ll are the main equations of the variable phase method. 
The method of variable phase has been developed by V.V. Babikov and others in the middle 


of 1960’s 


(see also references therein). Each of the two equations Eqs. fHHll - fHbll can 


be solved to high accuracy by using various numerical methods. In general, if we know the 
phase function di{r) for each r, then it is possible to determine the amplitude A£{r) of the 
wave function and the wave function \kfc(r) = AAr) exp{t6£{r)) itself. The corresponding 


equation for the amplitude ^^(r) takes the form 2 


dAi{r) 

dr 


Ai(r)V(r; a,b) r. ^ . /, \ . c- / \1 r ■ /r \ • e- / n 

— Jiykr) cos Oi[r) — ne[kr) sm oe[r) je[kr) sm oe[r) 


k 


+ni{kr) cos(5£(r) 


(47) 


Note that the variable phase method works very well for central potentials which rapidly 
vanish at r —)■ +cxo, e.g., for | V{r) |~ r"”, where n > 3 at r ^ +cx). This means that we 
can apply this method to determine the accurate wave functions of the continuous spectra 
in the held of the polarization potential, Eq. flddh . created by the neutral H atom which is 
assumed to be inhnitely heavy. It appears that the variable phase method provides a large 
number of advantages in applications to the problems related with the absorption/emission 
of radiation from the states of continuous spectra. 
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In the case of inverse bremsstrahlung the energy conservation law is written in the form 


nk^ 


+ U} = 


hk) 


, or hu = — -^ (48) 

2me 2me 2me 2me 

where fc* and kj are the wave numbers dehned as fca =| |= and pa is the momentum of 

the particle a and is the corresponding wave vector. In atomic units k^ = and ka = Pa, 
while Eq. fHHl) takes the form cu = \{p]—pi)- Also, from Eg. (1451) one hnds kf = ^kf + 

By using the variable phase method described above we can obtain the initial and hnal 
state wave functions, i.e. and 4/fc^(r), respectively. In the variable phase method 

these wave functions are represented as the products of their amplitudes and phase parts, 
i.e. exp(z5(>;fc,(r)) and exp{igi.kf(,r)). Let us assume 

that by using the procedure described above we have determined the incident and hnal wave 
functions. Now, these two wave functions can be used in the following computations of the 
cross-section of inverse bremsstrahlung. 

The explicit formula for the cross-section of inverse bremsstrahlung is (in relativistic 
units) 


d(T]ipjy 


ue 


e ■ pi^f 1^ dokd^pf 


(49) 


{2nYmePi 

where p* =| p* | is the ‘absolute’ momentum of the incident electron, p/ is the momentum 
of the hnal electron and Pi^f is the following matrix element 


Pi^f = -i 




(50) 


where and 4/^^ are the incident and hnal wave functions. By integrating Eg. (149 p over 
all directions of the hnal momentum one hnds the following formula for the cross-section (in 
atomic units): 


d(T\^p p — O: Ctf 


3 2^/ - 


e pi^f p doi^pjdpf 


° {2n)^Pi 

After a few steps of straightforward transformations one hnds from this formula 


2a^al {p} - pi)p) roo 

- - - I I iP! 


(61) 


(52) 


By using Egs. 
function: 
d'^i _ dAe-kiir) 


and (1471) we obtain the following expression for the derivative of the 4/ 

dSe-,ki{r) 


exp{t6e-ki) + ^At,ki{r) exp{t6i-ki)- 


dr dr 

I/(r;a,6)A£.fc(r)exp(z5<);fcJ 


dr 


(63) 


k 


ji{kr) cos Si{r) — ni{kr) sin 6i{r) je{kr) sin 6i{r) 


+ ne{kr) cos 6i{r) 


2 V ^Qj r 4 2 

- ^ ’ A(,.ki{r) exp{i6t,ki) \ji>{kr) cosSi(r) - ni(kr) sinSe(r 
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or, in other words, 


d'^i 

dr 


V{r;a,b) 

k 


je{kr) cos Si{r) — ni{kr) sin 6e{r) ji{kr) sin 6i{r) + ni{kr) cos dii^r) 


2V{r;a,b) r n \ ^ 

---M/j ji[kr) cosdi{r) — ni{kr) smo^^r 


The last equation does not include explicitly the amplitude function Note also that 

the last equation can be re-written in a number of different forms which are more convenient 
in some cases. 

Another simplihcation follows from the fact that the wavelength A of the light quantum 
which produces photodetachment of the H~ ion is > 12,000 A, while the effective ‘scatter¬ 
ing zone’ Rs for the process, Eq.(l2]) is less than 15 - 25 A. Therefore, we can apply the 
approximation based on small values of the dimensionless parameter ^ <C 1. This approx¬ 
imation allows one to replace the r—dependent phase factor 6^-kiif') in the wave function by 
its asymptotic value Si-kA+oo) determined at the inhnity. In other words, the phase of the 
wave function does not depend upon r, i.e. it is a constant. The wave function T£;fc.(r) is 
represented in the form = Ai-kiir) exp{z6£-ki), where dt = 5^{r = -|-cx)) is a constant, 

or almost a constant. This drastically simplihes numerical computations of all integrals 
needed in Eq. flS^ . Finally, it should be noted that based on the approach developed here 
we can evaluate the total cross-section of inverse bremsstrahlung for electrons scattered by 
the neutral hydrogen atoms in stellar photospheres. 


IX. CONCLUSION 

We have considered the absorption of infrared and visible radiation by the negatively 
charged hydrogen ions H~. Such an absorption of radiation plays an extremely important 
role in Solar and stellar astrophysics. In general, absorption of infrared and visible radiation 
by the H“ ions is mainly related with the photodetachment of these ions. The absorp¬ 
tion of radiation by electrons scattered in the held of the neutral hydrogen atom(s) also 
contributes to the radiation absorbtion coefficient at very large wavelength A > 12, OOOA 
(infrared radiation). In this study we discuss all details important in calculations of the 
photodetachment cross-section with the use of highly accurate (or truly correlated) wave 
functions of the two-electron (or three-body) H~ ion. Photodetachment cross-sections of 
the °°II~, ^11“ (or T“), (or D~) and ions are determined to high accuracy and for 
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a large number of photo-electron moment a/energies. Some additional details of our methods 


and computational results will be discussed in our next study 


29|. 


Right now, we can only note that our current calculations performed for the photode¬ 
tachment cross-section for the ion indicate a good agreement with the results obtained 
for this ion by other methods (see, e.g., [l^ - il6|). For the negatively charged hydrogen ions 
with the hnite nuclear masses calculations the photodetachment cross-sections with highly 
accurate wave functions have never been performed. In this sense our results for the 
(or T“), (or D“) and ions from Tables VI and VII are unique. Furthermore, this 
paper opens a new avenue for accurate computations of the photodetachment cross-section 
of all negatively charged hydrogen ions with the use of highly accurate wave functions. In 
future calculations there is no need to reduce the overall accuracy of the truly correlated 
wave functions by using different approximate procedures, e.g., by representing these wave 
functions in the form of natural orbital expansion(s). Currently, we are working on the prob¬ 
lem of incorporating extremely accurate wave functions with 3500 - 4000 basis functions (see 
Table I) in calculations of the photodetachment cross-sections. 

Appendix 


In this Appendix we test our variational wave functions constructed for the negatively 
charged hydrogen ions (°°II~, ^H~, ^H~ and in actual computations of the lowest order 
QED corrections for these two-electron systems. The lowest order QED corrections for the 
ground 1^5—state in the negatively charged hydrogen ion is determined by the formula 
(in atomic units) 


= -Qa 


^^—-2\na 


-(- Oi 


3 ^ L30 

rl64 14 


IniCo 
10 


(<^(reiv)) 


+ — \na- —S{S + 1) (5(ree)) - —a'’(—) 


7 


(54) 


L 15 3 3 ' Ott 

where a is the hne structure constant, Q is the nuclear charge (in atomic units) and S is the 
total electron spin. For the ground states in all H“ ions considered in this study we have 
Q = 1 and S' = 0 (singlet states) in Eq. fl54p . The notation IniFo in this formula stands for 
the Bethe logarithm 20|. To simplify numerical calculations in this Appendix one can take 
the value of In iFo from Table 3 in |5|. 

For the two-electron ions with hnite nuclear mass we need to evaluate the corresponding 
recoil correction to the lowest-order QED correction Eq. fIST)) . Such a correction is 
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calculated with the use of the following formula (in atomic units) 


- 


/ 2 


+ 


VM M + 1 


- Ilia - 4111 Ao 


3M L 3 


(<^(re7v)) 


7a3 1 

37rM 


(55) 


where M S> rUe is the nuclear mass. All expectation values in this equation must be 
determined for the real two-electron ions, i.e. ions with the hnite nuclear masses. The 
inverse mass ^ is a small dimensionless parameter which for the ions considered in this 
study is smaller than < 5 ■ 10“^. The last terms in both Eq.([5T) and Eq. fl5^ are usually 
called the Araki-Sucher terms, or Araki-Sucher corrections 30|, 3]|, since this correction was 
obtained and investigated for the hrst time in papers by Araki and Sucher. The expectation 
value of the term (-|-) is singular, i.e. it contains the regular and non-zero divergent parts. 

^ee 

General theory of the singular three-body integrals was developed in our earlier works (see, 
e.g., B and references therein). In particular, in we have shown that the (^) expectation 


value is determined by the formula 


+ 47r((5(rgg)) (56) 

^ee ^ ee 

where (-|-)r is the regular part of this expectation value and ((5(ree)) is the expectation value 

^ee 

of the electron-electron delta-function. Briefly, we can say that the overall contribution of 
the signular part of the ^ operator is reduced to the expectation value of the corresponding 

^ee 

delta-function. 

It follows from the formulas Eqs. (l54p and (l55ll that highly accurate evaluation of the 
lowest-order QED correction is a very effective test of the non-relativistic wave functions of 
the negatively charged hydrogen ions. After this test the bound state wave functions can 
be used in numerical calculations of the photodetachment cross-section(s) and absorption 
coefficients of the negatively charged hydrogen ions. Indeed, the formulas Eqs. dSTj) and 
(l5^ contain both the electron-nuclear and electron-electron delta-functions. To prove the 
correctness of these values and approximately evaluate the total number of stable decimal 
digits in each of the delta-functions one need to compute the corresponding cusp values and 
compare them with the expected (or predicted) values. This essentially coincides with the 
standard test for the non-relativistic wave functions. In addition to these delta-functions 
Eqs. fl54p and fl5^ include the (^) expectation value which is singular. Accurate numerical 
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evaluation of the (-3—} and (-|-) expectation values is not an easy task (see, e.g., |5[), but 

^eJV 

it provides another effective test for the non-relativistic wave function. Our results for the 
lowest-order QED correction for the and ions can be found in Table V 

(in atomic units). The lowest-order QED corrections and AE^^ determined for 

each of these ions with the hnite nuclear masses are the most accurate values to-date. In 
reality, they coincide in seven hrst decimal digits with the corresponding values from j^, 
where we used less accurate variational wave functions. Very likely, the lowest-order QED 
corrections are of interest by themselves, but we have calculated them here only to prove 
very high quality of our non-relativistic wave functions. 
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TABLE I: The non-relativistic total energies E (in a.u.) determined with different number of 
basis functions N, the asymptotic values of the total energies E{N = oo) (in a.u.) and ioniza¬ 
tion/detachment potential xi (iii cE) of the hydrogen negatively charged ions (isotopes) in atomic 
units. 


K 

E(°°H-) 


e(3h-) 

3500 

-0.527751 016544 377196 590213 


-0.527649 048201 920733 538766 

3700 

-0.527751 016544 377196 590333 


-0.527649 048201 920733 538885 

3840 

-0.527751 016544 377196 590389 


-0.527649 048201 920733 538983 

4000 

-0.527751 016544 377196 590446 


-0.527649 048201 920733 539050 

E{N = oo) 

-0.527751 016544 377196 59075(10) 


-0.527649 048201 920733 539355(10) 

Xi 

-0.755143 903366 701124 238 


-0.754843 900893 517142 964 

K 

E{^R-) 


E(iH-) 

3500 

-0.527598 324689 706528 595827 


-0.527445 881119 767477 071142 

3700 

-0.527598 324689 706528 595960 


-0.527445 881119 767477 071261 

3840 

-0.527598 324689 706528 596044 


-0.527445 881119 767477 071359 

4000 

-0.527598 324689 706528 596110 


-0.527445 881119 767477 071425 

E{N = oo) 

-0.527598 324689 706528 596355(10) 


-0.527445 881119 767477 071665(10) 

Xi 

-0.754694 721951 430623 253 


-0.754246 603605 793792 368 
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TABLE II: Convergence of the expectation values of some electron-nuclear and electron-electron 
properties for the ion. All values are shown in atomic units and N is the total number of 

basis functions used. 



{reN) 

(rl) 

(^eTv) 

3500 

2.710178278444420365208 

25.20202529124033189252 

1.11666282452545228 

3700 

2.710178278444420365265 

25.20202529124033189563 

1.11666282452544953 

3840 

2.710178278444420365286 

25.20202529124033189660 

1.11666282452544421 

4000 

2.710178278444420365301 

25.20202529124033189738 

1.11666282452543860 


T7V 


(reAT • Fee) 

3500 

0.6451445424122193894344 

0.26387550827218859829282 

12.60101264562016594726 

3700 

0.6498715811920881669355 

0.26387550827218859829304 

12.60101264562016594781 

3840 

0.6498715811920881669351 

0.26387550827218859829337 

12.60101264562016594830 

4000 

0.6498715811920881669349 

0.26387550827218859829358 

12.60101264562016594869 


(“^The notation TsN stands for the (cos(r 3 i^r 2 i)) = (cos(reAr^ree)) expectation value. Note also, 
that for all two-electron atomic systems 2TeN + Tee = 1 + 4(/), where (/) is the expectation value 

mentioned in Table III. 
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TABLE III: The expectation values of some propeties (in atomic units) for the °°H ion. 


(^eiv) 

{r-J) 

(^eTv) 


1.11666282452542572 

0.15510415256242466 

0.6832617676515272224 

0.311021502214300052 

(reN) 

(^ee) 

{r^N) 

(rl) 

2.7101782784444203653 

4.4126944979917277211 

11.913699678051262274 

25.202025291240331897 

(^eAr) 

{rle) 

(rtN) 

(rte) 

76.02309704902717911 

180.60560023017477483 

645.144542412219375 

1590.09460393948530 

(N2r'3i]"^) 

{[reNTeeV^) 

([r32r3ir'2i]“^) 

((5(reeAr)) 

0.38262789034020545 

0.25307756706456687 

0.20082343962918944 

5.129778775490-10-3 

(<^(reiv)) 

(a) 

^eN 

(<^(reiv)) 

(a) 

Ree 

0.1645528728473590 

-1.00000000001778 

2.737992126104611-10-3 

0.500000002446 

TeN 

'^ee 

if) 

(r3i • r32/rii) 

0.6498715811920881669 

-0.1051476935659779011 

0.048648867204549608205 

-0.4642618530806317014 



(Ve • Ve) 

(Ve • Vat) 

0.2638755082721885983 

0.560630798396681918 

0.0328797818523047217 

-5.60630798396681918 

{VeN ■ Tee) 

{irJ)R) 

{iree)R) 

(’’eTv) 

12.601012645620165948 

-3.43559485054432 

0.064307887285283 

-1.36776246468689 


(a)The expected cusp values (in a.u.) for the ion are = —1.0 and t'ee = 0.5 (exactly). 
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TABLE IV: Convergence of the expectation values of the electron-nuclear and electron-electron 
delta-functions {{5{re]\f)) and ((5(ree))) and regular parts od the inverse cubic expectation values 
and {r~J^)R) for the and ions. All values are shown in atomic units and N 

is the total number of basis functions used. 


■ 

(tritium) 

2H (deuterium) 

(protium) 

■ 

(<^(reiv)) 

(<^(re7v)) 

(<^(re7v)) 

3500 

0.16446163681098 

0.16441626335994 

0.16427994412837 

3700 

0.16446163681128 

0.16441626336024 

0.16427994412867 

3840 

0.16446163681069 

0.16441626335966 

0.16427994412847 

4000 

0.16446163681105 

0.16441626336002 

0.16427994412844 

■ 

(<^(ree)) 

(<^(ree)) 

(<^(ree)) 

3500 

2.735845627669-10-3 

2.734778360513-10-3 

2.731572792726-10-3 

3700 

2.735845627652-10-3 

2.734778360496-10-3 

2.731572792710-10-3 

3840 

2.735845627712-10-3 

2.734778360555-10-3 

2.731572792769-10-3 

4000 

2.735845627776-10-3 

2.734778360619-10-3 

2.731572792833-10-3 

■ 




3500 

-3.43332137343396 

-3.43219080091021119 

-3.4287944147910614 

3700 

-3.43332137346264 

-3.43219080093883454 

-3.4287944148194641 

3840 

-3.43332137341561 

-3.43219080089185737 

-3.4287944147726374 

4000 

-3.4333213734460 

-3.43219080092215666 

-3.4287944148027678 

■ 

{rIe)R 

{ree)R 

{ree)R 

3500 

6.427894892817-10-2 

6.426455565395-10-2 

6.422130608689-10-2 

3700 

6.427894893018-10-2 

6.426455565595-10-2 

6.422130608885-10-2 

3840 

6.427894892413-10-2 

6.426455564991-10-2 

6.422130608281-10-2 

4000 

6.427894892942-10-2 

6.426455564520-10-2 

6.422130608203-10-2 
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TABLE V: The Bethe logarithm (in a.u.) and lowest order QED corrections and (in 

MHz) for the ground 1^5—state(s) in the and ions. 



OOH- 

(tritium) 

(deuterium) 

(protium) 

IniLo {a.u.) 

2.993004415 

2.993011414 

2.993014897 

2.993025369 


8215.203465 

8210.661159 

8208.402174 

8201.615291 

AE^^^ 

8215.203465 

8207.191206 

8203.207102 

8191.239514 
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TABLE VI: Photodetachment cross-section (in cm?) of the negatively charged hydrogen ion in 
their ground states. N is the total number of basis function in the trial wave function and 
E is the total energy (in a.u.) for this wave function. The final hydrogen atom is formed in the 
ground Is-state. The notation Eg stands for the kinetic energy of the emitted photo-electron (in 

a.u.). 



oog- 

oog- 

3h- 

2h- 

iR- 

Eh- 

-0.5277751016544153 

-0.5277751016544252 

-0.527649048201671 

-0.527598324689580 

-0.527445881119594 

Ee 

a{Ee){N = 300) 

<j{E^){N = 350) 

cr{Eg){N = 350) 

(T{Ef,){N = 350) 

a{Ee)(N = 350) 

0.0000125 

0.2959861281153E-17 

0.2959854603169E-17 

0.2959608363575E-17 

0.2959527360520E-17 

0.2959202302718E-17 

0.0001125 

0.8795361962978E-17 

0.8795347268968E-17 

0.8794584542552E-17 

0.8794308561005E-17 

0.8793268085947E-17 

0.0003125 

0.1438357679058E-16 

0.1438356148152E-16 

0.1438220490456E-16 

0.1438164996520E-16 

0.1437971473501E-16 

0.0006125 

0.1957704797715E-16 

0.1957703477352E-16 

0.1957494491642E-16 

0.1957401257879E-16 

0.1957094252745E-16 

0.0010125 

0.2425248779858E-16 

0.2425247518253E-16 

0.2424946573180E-16 

0.2424806431545E-16 

0.2424359742878E-16 

0.0015125 

0.2831673151310E-16 

0.2831671287514E-16 

0.2831258986430E-16 

0.2831063012238E-16 

0.2830450229469E-16 

0.0021125 

0.3170965737411E-16 

0.3170962354548E-16 

0.3170421932053E-16 

0.3170161206699E-16 

0.3169358006790E-16 

0.0028125 

0.3440414647798E-16 

0.3440409146283E-16 

0.3439728216636E-16 

0.3439394843234E-16 

0.3438381596859E-16 

0.0036125 

0.3640350418486E-16 

0.3640342883386E-16 

0.3639514145850E-16 

0.3639102647578E-16 

0.3637866660179E-16 

0.0045125 

0.3773699512169E-16 

0.3773690666895E-16 

0.3772712159963E-16 

0.3772220418799E-16 

0.3770757332968E-16 

0.0055125 

0.3845424767129E-16 

0.3845415662293E-16 

0.3844290725784E-16 

0.3843720233167E-16 

0.3842034419486E-16 

0.0066125 

0.3861926168889E-16 

0.3861917835696E-16 

0.3860654758202E-16 

0.3860010314572E-16 

0.3858114331217E-16 

0.0078125 

0.3830463436882E-16 

0.3830456649096E-16 

0.3829068002096E-16 

0.3828357090908E-16 

0.3826270448851E-16 

0.0091125 

0.3758644736645E-16 

0.3758639917081E-16 

0.3757141647769E-16 

0.3756373675088E-16 

0.3754121197964E-16 

0.0105125 

0.3654007536656E-16 

0.3654004777150E-16 

0.3652415173423E-16 

0.3651600703972E-16 

0.3649210748110E-16 

0.0120125 

0.3523701334398E-16 

0.3523700477057E-16 

0.3522039120618E-16 

0.3521189189400E-16 

0.3518691928177E-16 

0.0136125 

0.3374269418098E-16 

0.3374270152781E-16 

0.3372556954001E-16 

0.3371682494788E-16 

0.3369108401747E-16 

0.0153125 

0.3211518525441E-16 

0.3211520475640E-16 

0.3209774852243E-16 

0.3208886257917E-16 

0.3206264861165E-16 

0.0171125 

0.3040460861272E-16 

0.3040463646812E-16 

0.3038703877562E-16 

0.3037810689795E-16 

0.3035169622355E-16 

0.0190125 

0.2865311628737E-16 

0.2865314904517E-16 

0.2863557664068E-16 

0.2862668378989E-16 

0.2860032716430E-16 

0.0210125 

0.2689526056870E-16 

0.2689529532263E-16 

0.2687789595522E-16 

0.2686911566400E-16 

0.2684303422298E-16 

0.0231125 

0.2515861990078E-16 

0.2515865435015E-16 

0.2514155528590E-16 

0.2513294944699E-16 

0.2510733283533E-16 

0.0253125 

0.2346456740216E-16 

0.2346459983308E-16 

0.2344790756524E-16 

0.2343952678096E-16 

0.2341453301346E-16 

0.0276125 

0.2182909599163E-16 

0.2182912521773E-16 

0.2181292608860E-16 

0.2180481040352E-16 

0.2178056702623E-16 

0.0300125 

0.2026363869442E-16 

0.2026366398150E-16 

0.2024802545112E-16 

0.2024020531988E-16 

0.2021681148092E-16 

0.0325125 

0.1877584342631E-16 

0.1877586441441E-16 

0.1876083674625E-16 

0.1875333412069E-16 

0.1873086321621E-16 

0.0351125 

0.1737027795553E-16 

0.1737029458465E-16 

0.1735591275960E-16 

0.1734874220901E-16 

0.1732724483298E-16 

0.0378125 

0.1604905302160E-16 

0.1604906546358E-16 

0.1603535117581E-16 

0.1602852097992E-16 

0.1600802796753E-16 

0.0406125 

0.1481236028999E-16 

0.1481236888700E-16 

0.1479933251391E-16 

0.1479284568916E-16 

0.1477337109967E-16 

0.0435125 

0.1365892760874E-16 

0.1365893281871E-16 

0.S64657529017E-16 

0.1364043052627E-16 

0.1362197445888E-16 





















TABLE VII: Photodetachment cross-section (in cm^) of the negatively charged hydrogen ions in 
their ground 1^5—states for intermediate (kinetic) energies of photo-electron Eg. All notations are 
the same as in the previous Table. 



oog- 

oog- 

3h- 

2H- 

iR- 

Ee 

<j{Ee){N = 300) 

<j{E^){N = 350) 

cr{Ee){N = 350) 

(T{Ef,)(N = 350) 

a(Ef,){N = 350) 

0.0465125 

0.1258639757846E-16 

0.1258639992733E-16 

0.1257471444683E-16 

0.1256890694264E-16 

0.1255145811860E-16 

0.0496125 

0.1159163734473E-16 

0.1159163738578E-16 

0.1158061097734E-16 

0.1157513317341E-16 

0.1155867124808E-16 

0.0528125 

0.1067098826764E-16 

0.1067098654801E-16 

0.1066060141832E-16 

0.1065544362939E-16 

0.1063994122132E-16 

0.0561125 

0.9820464104618E-17 

0.9820461139031E-17 

0.9810695855883E-17 

0.9805846809755E-17 

0.9791271250301E-17 

0.0595125 

0.9035905851910E-17 

0.9035902103435E-17 

0.9026732599035E-17 

0.9022179890367E-17 

0.9008494708252E-17 

0.0630125 

0.8313100637491E-17 

0.8313096504894E-17 

0.8304496917338E-17 

0.8300227385501E-17 

0.8287393547857E-17 

0.0666125 

0.7647871188006E-17 

0.7647866999040E-17 

0.7639810367876E-17 

0.7635810408153E-17 

0.7623787342894E-17 

0.0703125 

0.7036141501167E-17 

0.7036137510947E-17 

0.7028596343733E-17 

0.7024852160681E-17 

0.7013598586415E-17 

0.0741125 

0.6473983501310E-17 

0.6473979894812E-17 

0.6466926621475E-17 

0.6463424432608E-17 

0.6452899021906E-17 

0.0780125 

0.5957648673312E-17 

0.5957645571010E-17 

0.5951052902029E-17 

0.5947779099168E-17 

0.5937940999571E-17 

0.0820125 

0.5483587974654E-17 

0.5483585440313E-17 

0.5477426643500E-17 

0.5474367918339E-17 

0.5465177153440E-17 

0.0861125 

0.5048462721451E-17 

0.5048460771207E-17 

0.5042709880386E-17 

0.5039853319997E-17 

0.5031271090582E-17 

0.0903125 

0.4649148629125E-17 

0.4649147241080E-17 

0.4643779208428E-17 

0.4641112365692E-17 

0.4633101268488E-17 

0.0946125 

0.4282734755358E-17 

0.4282733879064E-17 

0.4277724681340E-17 

0.4275235624615E-17 

0.4267759803810E-17 

0.0990125 

0.3946518733255E-17 

0.3946518298616E-17 

0.3941845006054E-17 

0.3939522351764E-17 

0.3932547596255E-17 

0.1035125 

0.3637999386976E-17 

0.3637999312224E-17 

0.3633640127535E-17 

0.3631473058861E-17 

0.3624966856151E-17 

0.1081125 

0.3354867581199E-17 

0.3354867779736E-17 

0.3350802053473E-17 

0.3348780327823E-17 

0.3342711883700E-17 

0.1128125 

0.3094995961090E-17 

0.3094996347097E-17 

0.3091204574588E-17 

0.3089318522133E-17 

0.3083658753701E-17 

0.1176125 

0.2856428083343E-17 

0.2856428576147E-17 

0.2852892380277E-17 

0.2851132895481E-17 

0.2845854403773E-17 

0.1225125 

0.2637367314499E-17 

0.2637367841750E-17 

0.2634069945105E-17 

0.2632428472874E-17 

0.2627505500110E-17 

0.1275125 

0.2436165774447E-17 

0.2436166274218E-17 

0.2433090463444E-17 

0.2431558981225E-17 

0.2426967356852E-17 

0.1326125 

0.2251313525848E-17 

0.2251313947789E-17 

0.2248445032448E-17 

0.2247016029429E-17 

0.2242733108233E-17 

0.1378125 

0.2081428150094E-17 

0.2081428455791E-17 

0.2078752223555E-17 

0.2077418677866E-17 

0.2073423272820E-17 

0.1431125 

0.1925244803977E-17 

0.1925244966678E-17 

0.1922748136311E-17 

0.1921503490986E-17 

0.1917775802922E-17 

0.1485125 

0.1781606815719E-17 

0.1781606819574E-17 

0.1779276992832E-17 

0.1778115131145E-17 

0.1774636676911E-17 

0.1540125 

0.1649456852113E-17 

0.1649456691077E-17 

0.1647282304415E-17 

0.1646197525071E-17 

0.1642951065436E-17 

0.1596125 

0.1527828668478E-17 

0.1527828344977E-17 

0.1525798621745E-17 

0.1524785614298E-17 

0.1521755082618E-17 

0.1653125 

0.1415839438357E-17 

0.1415838961857E-17 

0.1413943865477E-17 

0.1412997686281E-17 

0.1410168118652E-17 

0.1711125 

0.1312682649388E-17 

0.1312682034922E-17 

0.1310912223441E-17 

0.1310028272342E-17 

0.1307385739798E-17 

0.1770125 

0.1217621544363E-17 

0.1217620811106E-17 

0.1215967593385E-17 

0.1215141591308E-17 

0.1212673134477E-17 
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TABLE VIII: Photodetachment cross-section (in cm?) of the negatively charged hydrogen ions in 
their ground I ^ S'—states for kinetic energies of photo-electron which are far from the maximum 
of the cross-section. All notations are the same as in the two previous Tables. 



oog- 

oog- 

3h- 

2h- 

iR- 

Ee 

a{Ee){N = 300) 

a{Ee){N = 350) 

a{Ee){N = 350) 

cr{Ee){N = 350) 

a{Ee)(N = 350) 

0.1830125 

0.1129983081672E-17 

0.1129982251607E-17 

0.1128437545329E-17 

0.1127665512883E-17 

0.1125359079348E-17 

0.1891125 

0.1049152386308E-17 

0.1049151483018E-17 

0.1047707774635E-17 

0.1046986011787E-17 

0.1044830396360E-17 

0.1953125 

0.9745676610310E-18 

0.9745667086430E-18 

0.9732170153305E-18 

0.9725420821878E-18 

0.9705268701701E-18 

0.2016125 

0.9057155267702E-18 

0.9057145490617E-18 

0.9044523827004E-18 

0.9038210814197E-18 

0.9019365949053E-18 

0.2080125 

0.8421267615897E-18 

0.8421257812580E-18 

0.8409451144388E-18 

0.8403544722382E-18 

0.8385917194393E-18 

0.2145125 

0.7833724083041E-18 

0.7833714463962E-18 

0.7822666804163E-18 

0.7817139336736E-18 

0.7800645612302E-18 

0.2211125 

0.7290602220173E-18 

0.7290592975088E-18 

0.7280252323105E-18 

0.7275078117097E-18 

0.7259640599065E-18 

0.2278125 

0.6788314302517E-18 

0.6788305597618E-18 

0.6778623657516E-18 

0.6773778824334E-18 

0.6759325432339E-18 

0.2346125 

0.6323577798990E-18 

0.6323569775261E-18 

0.6314501692030E-18 

0.6309964018707E-18 

0.6296427796583E-18 

0.2415125 

0.5893388468558E-18 

0.5893381241154E-18 

0.5884885354307E-18 

0.5880634183601E-18 

0.5867952932704E-18 

0.2485125 

0.5494995858601E-18 

0.5494989517092E-18 

0.5487027130754E-18 

0.5483043249749E-18 

0.5471159186943E-18 

0.2556125 

0.5125880996827E-18 

0.5125875606145E-18 

0.5118410774791E-18 

0.5114676311479E-18 

0.5103535750516E-18 

0.2628125 

0.4783736084176E-18 

0.4783731686084E-18 

0.4776731015090E-18 

0.4773229342409E-18 

0.4762782397414E-18 

0.2701125 

0.4466446011363E-18 

0.4466442626341E-18 

0.4459875086486E-18 

0.4456590733307E-18 

0.4446791043100E-18 

0.2775125 

0.4172071536235E-18 

0.4172069165616E-18 

0.4165905920736E-18 

0.4162824489309E-18 

0.4153628961306E-18 

0.2850125 

0.3898833972676E-18 

0.3898832600957E-18 

0.3893046847916E-18 

0.3890154937331E-18 

0.3881523509801E-18 

0.2926125 

0.3645101254581E-18 

0.3645100851790E-18 

0.3639667672004E-18 

0.3636952807213E-18 

0.3628848228715E-18 

0.3003125 

0.3409375250224E-18 

0.3409375774283E-18 

0.3404271996037E-18 

0.3401722562045E-18 

0.3394110186911E-18 

0.3081125 

0.3190280213342E-18 

0.3190281612350E-18 

0.3185485683200E-18 

0.3183090864030E-18 

0.3175938462811E-18 

0.3160125 

0.2986552267348E-18 

0.2986554481751E-18 

0.2982046350311E-18 

0.2979796072405E-18 

0.2973073656244E-18 

0.3240125 

0.2797029828421E-18 

0.2797032793025E-18 

0.2792793799499E-18 

0.2790678679205E-18 

0.2784358337181E-18 

0.3321125 

0.2620644881721E-18 

0.2620648527520E-18 

0.2616661302361E-18 

0.2614672597225E-18 

0.2608728345722E-18 

0.3403125 

0.2456415032806E-18 

0.2456419288630E-18 

0.2452667658745E-18 

0.2450797222305E-18 

0.2445204865560E-18 

0.3486125 

0.2303436263476E-18 

0.2303441057424E-18 

0.2299909959384E-18 

0.2298150199219E-18 

0.2292887200228E-18 

0.3570125 

0.2160876327733E-18 

0.2160881588421E-18 

0.2157556988156E-18 

0.2155900826944E-18 

0.2150946187971E-18 


38 





















